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Abstract 

Using total variation based energy minimisation we address the recovery of a blurred 
(convoluted) one dimensional (ID) bar code. We consider functionals defined over all 
possible bar codes with fidelity to a convoluted signal of a bar code, and regularised 
by total variation. Our fidelity terms consist of the distance either directly to the 
measured signal or preceded by deconvolution. Key length scales and parameters are 
the X-dimension of the underlying bar code, the size of the supports of the convolution 
and deconvolution kernels, and the fidelity parameter. For all functionals, we establish 
parameter regimes (sufficient conditions) wherein the underlying bar code is the unique 
minimiser. We also present some numerical experiments suggesting that these sufficient 
conditions are not optimal and the energy methods are quite robust for significant 
blurring. 
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1 Introduction and notation 

A one- dimensional bar code is a finite series of alternating black bars and white spaces 
with varying widths (Figure [T]). The so-called X -dimension of the bar code is the width 
of the narrowest bar or space. A bar code scanner can use, for example, light detectors or 
photoreceptors similar to those used by a digital camera to pick up light and dark areas 
of the bar code, and produce a continuous signal associated with the darkness distribution 
across the span of the bar code (c.f. |10)). One is left with an approximation to the bar 
code which depending on the scanner and the way in which the scan was taken (distance, 
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Figure 1: Left: A standard ID 12-digit Universal Product Code (UPC) bar code. The 
12 numerical digits readable to the human eye are encoded in the bar spacings. Right: 
A bar code scanner in which the black bars absorb light, while the white bars reflect it. 
Photo diodes turn the reflected light into an electrical signal, which may be both blurred 
and noisy. This signal is then converted into digital pulses. 

vibrations etc.) can be blurred and noisy, (see Figure [T]). Thus in this article we deal 
with the general question: Given a blurred and possibly noisy signal f associated with a 
bar code, how can we deblur and denoise effectively to reconstruct the bar code? 

Standard commercial deblurring techniques are often based upon edge detection, for 
example, finding local extrema of /' which hopefully correspond to the discontinuities (i.e. 
interfaces) of the original bar code. As noted in |i4j, this presents several difficulties: (i) 
the process is highly unstable to small changes in the signal, for example to the presence of 
noise; (ii) points associated with local extrema of /' are only crude approximations to the 
true locations of the bars with convolution tending to distort these points if the standard 
deviation of the convolution kernel — or if the kernel is compactly supported the size of its 
support — is comparable to the X-dimension of the bar code; (iii) if the standard deviation 
or support size of the kernel is very large compared to the X-dimension, some edges in the 
bar code may have no corresponding extrema of /' at all. 

In this article, we take a different approach based upon energy minimisation involving a 
total variation (TV) regularisation - a method introduced by Rudin, Osher and Fatemi in 

. One of the many advantages of this approach is that energy minimisation is a process 
stable with respect to small changes in the input signal. 

In the context of bar code reconstruction, this type of energy minimisation was first 
proposed by Fadil Santosa, and analysed by Selim Esedoglu [1] (see also [H]). Here we take 
a similar approach but with an important difference: We are interested in directly testing 
the merits of energy minimisation for TV-based functionals by considering ansdtze for / 
which involve convolutions of a bar code with certain blurring kernels. The functionals, 
blurring kernels, and admissible classes of bar codes possess certain length scale parameters. 
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In [1] Esedoglu shows existence of solutions for a variety of functionals and then proceeds to 
numericahy test these functionals via an algorithm which approximates the actual length 
scale parameter of the blurring kernel and uses this information to reconstruct the clean bar 
code signal. Here we fix these length scales as parameters and ask under what conditions 
can we insure that energy minimisation gives back the underlying bar code. 

We begin by introducing some notation. A bar code is given by a function u of bounded 
variation (c./. [6]) with suppu C [0, 1] taking on the values and 1 a.e., i.e. we consider a 
subset of the space BV{M; {0, 1}). In particular the general admissible set for bar codes is 

S := |u G BV{R; {0, 1}) : u = a.e. in [0, 1]' 

Modulo a set of measure zero, the set {x G [0, 1] : u{x) = 1} consists of a finite number 
of disjoint non-empty intervals called bars. We denote these non-empty bar intervals by 
bi with length \bi\. Similarly the intervals of nonempty white spaces (i.e. intervals in [0, 1] 
where u = 0) are denoted by Wi. In this paper, signals will always be generated by bar 
codes whose X-dimension is bounded below. That is, we assume there exists a constant 
00 > such that the minimum width of these bars and spaces is a priori bounded below 
by UJ. Our space of generating bar codes is thus 

Buj := < u £ B : Vi min{|6j|, \wi\} > uj 



In some calculations it is useful, though harmless, to assume that we know whether the 
bar code starts (and/or ends) with a bar or a space. To this end, we define for i,j G {0, 1}, 
the sets 

B^^ := |u G ;B : there exists xi, X2 such that u = i on [0, xi] and u = j on [x2, 

Bljj := |u G : there exists xi, X2 such that u = i on [0, xi] and u = j on [x2, 

As a useful addition to our terminology we will call the transition from a bar to a space or 
vice versa an interface. This means that /jg is equal to the number of interfaces in u 
(where interfaces at x = and x = 1 are included). The notation xs will be used for the 
characteristic function of a given set S. 

We approach the blurring via convolution on a length scale a > (c.f. [71[8lll2lH]) with 
a symmetric unimodal kernel of unit mass. Different results hold for kernels of varying 
generality and we will discuss these kernels shortly. For the moment, let the kernel (pa 
denote a symmetric probability distribution on M with size (for example, its standard 
deviation or if compactly supported, half the size of its support) a. Given a bar code 
zeB we will assume observed signals of the form: 

POO 

fa{x) := {(i)a*z){x) = / (t)a{x-y)z{y)dy. 
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where o" > is fixed (note that 4>o = 6, the Dirac delta distribution, and thus fo = z). 
For u & B: 

• We consider a fidehty term which compares u directly with the observed signal: 



FiH := / \u'\ + X\\u - /, 

JR 
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• {deconvolution/dehlurring functional) Under the belief that the signal involves a con- 
volution with a known kernel, we may incorporate this convolution into the structure 
of the fidelity term and consider: 

F2{u) := / \u \ + X\\(j)a*U - fa\\l2(^). 

If is known, one may ask as to the merits of energy minimisation as we could simply 
Fourier transform the observed signal and divide by the Fourier transform of (pa 
to recover the Fourier transform of z. However, this process is highly unstable with 
respect to small perturbations and in practice, there is always some noise associated 
with the observed signal f^. Directly solving back for z is analogous to solving the 
heat equation backwards. Energy minimisation provides a stable numerical approach 
to deblurring and indeed denoising. 

• (blind deconvolution/dehlurring functional) Assuming no knowledge of a, we may 
deconvolute with a kernel of similar type but with size p: 

F^iu) := / + X\\(l)p*u - faWhm- 

Note that -Fi and F2 are special cases of F3 with p equal to and a respectively, and that in 
our notation, the dependence of Fi on A, z, (pa and (pp is suppressed. It is straightforward 



to see (cf. Lemma 2.1), that for all parameters and all generating bar codes z, a minimiser 
of Fi(u) over u €z B exists. In this article, we examine the following questions: For what 
values of the parameters \,uj,a and p, is the minimiser of Fi known, and in particular, 
when is it the underlying bar code z7 Our results consist of two parts: First off we present 
a simple result that if A is sufficiently small, the unique minimiser is simply u = 0. Part 2 
deals with sufficient conditions for when the unique minimiser is the underlying bar code 
z. 

Note that if A is sufficiently small, n = is the unique minimiser in B (i.e. an empty 
bar code). This is clearly the case if 



A<Ao :=2/||^"2 



L2' 



since any nontrivial bar code has a minimum total variation of 2. Sufficient and necessary 
conditions on A for li = to be the unique minimiser in BV{M.) are given in the first 
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parts of the Theorems |1.2| and Corollary 1.3 by exploiting a method adopted from [9j 
wherein the following dual norm associated with BV is used: 



sup 



fv 



: V G L^{R) and / \v'\ < 1 L (1) 



This threshold for a trivial minimiser is given by A < y , p > (parts 1 of the 

theorems below) and is lower than Aq- To see this note that for any / G -L^(M; [0, 1]) with 
compact support we may take I to be a bounded subset of M such that supp (j)p * f C I. 
Since v = ^x/ is an admissible function in the definition of * we have 

Therefore 

' <JffY'< ' 
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For A between this threshold and Aq, = is the unique minimiser in B but not in BV{M). 

The case where z is the unique minimiser is more subtle and depends critically on both 
the size and particular nature of the blurring kernel. We make two assumptions here: 

1. We restrict our attention to kernels with compact and small (with respect to the 
X-dimension w) support. 

2. We further consider unimodal, symmetric kernels in IC defined below and explicit 
regimes are computed using a prototype of such a kernel, the hat function defined by 

{(1 - x/a)/a if < X < fj, 

(1 + x/a)/a if - a < X < 0, 
if \x\ > a. 

The class /C is defined by 

/C := {(pa G L^(M) : 3(7 > (paix) = p{-x,a)x[-a,o](.x) + p{x,cr)x[o,a]{x) 

for a non-negative function p : [0, a] x (0, oo) monotonically decreasing 

in X, and / p{x, a) dx = - \ . (2) 



We will consistently use a subscript as in (j)^ to indicate the value of the parameter 
a in the definition of fC, i.e. (/>p G /C is in the subset of /C where a = p. 

Note in particular that if 4)^ G L-'^(M), we have (pa *u a L^(M) for all u £ B. 

It is also convenient to consider a subclass of fC 

/C3 := {(pa G /C : (/)o- £ Cc(J^),p continuously differentiable in a and ^ holds} , 
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where the non-obvious condition ^ is 

Vr G {0,a], Vc > 2a, Vx G [0,c] : 

r r d 

J{a,T,x,c) := J J —p{y,T)[(j)a{y-w) + (/)a{y + w)]dwdy<0. (3) 



This condition insures a certain monotonicity property (c.f. Lemma 2.13) of double 
convolutions with bar codes. As we show in Appendix |B] a sufficient condition for ([s]) 
to be satisfied is if for each r G (0, a] 

d 

(a) either ^p(x,r) is monotonically increasing in x and J'{a,T,0, c) < for all 



d 

(b) or —p{x, r) is monotonically decreasing in x and J' {a, r, |, c) < for all c > 2a. 



In the same appendix we show that the hat function (p^ satisfies condition (a) above 
for each r G (0,ct]. In practical applications using kernels that do not satisfy either 
(a) or (b), condition Q can be tested numerically. 

Assumption [l] is rather important and indeed restrictive as it limits the possible effect of 
blurring. For F2 and we require a, p < lu/2, insuring no interactions between neigh- 
bouring bars. For Fi, the condition is slightly less restrictive, namely a < u. The second 
assumption, particularly, the specification of the hat function is more for convenience. A 
crucial step in our direct and rather brute-force approach is to assume a minimiser with 
a certain structure and directly construct competitors which differ on a set bar or space. 
For this step, one can explicitly calculate a regime wherein such a competitor exists, and 
for simplicity we have performed the calculations for the hat function (which were greatly 



simplified by the use of Maple). We discuss modifications for other kernels in Remark 1.4 
below. Let us now state our results. 

Theorem 1.1. The following hold: 

1. Let (f)(j ^ K,. u = is the unique minimiser for Fi over BV{M.) iff \\fa\\* < 

2. Let (pa- = (pa o,nd z £ B^j- If a and A satisfy 

2 2 

a < CO and -a + — < uj, (4) 

3 A 

then u = z is the unique minimiser of Fi over B. 
Theorem 1.2. The following hold: 

1. Let (pfj £ K,. u = is the unique minimiser for F^ over BV{M.) iff \\cpp * fa\\* ^ jX' 
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2. Let 



and z G B^J for some i,j G {0, 1}. Let a < p < %. If X, p, and a satisfy 



2 1 

+ 



+ bpa"^ + 17 p^] < co, 



A ' 15p2 

then u = z is the unique minimiser of F3 over B^^ . 



(5) 



Note that the left hand side of ([5]) is increasing as a function of (real and positive) p 
and a. By taking p = a va. Theorem 1.2, we obtain 



Corollary 1.3. The following hold: 

1. Let (pa ^ K.- u = is the unique minimiser for F2 over BV{M) iff 

||</><7 * /all* < ^• 

2. Let (pa = (pa and z G for some i,j G {0, 1}. If a < ^ and A > satisfy 

2 21 

A + 15^<"' 
then u = z is the unique minimiser of F2 over B^^ . 



(6) 



Remark 1.4. Extensions to other kernels We remark on extensions of parts 2 of 
Theorem |1.1[ Theorem |1.2[ and Corollary 1.3 to more general kernels in /C. Their proofs 
consist of two steps: (i) First is to establish that any minimiser of Fi or F^ distinct from 
z must have strictly less interfaces than z (note that this is trivially satisfied for F2 since 
z uniquely minimises the fidelity term). For Fi, this is a simple consequence of vanishing 
first variation (Lemma 2.6), and holds for any kernel in fC. For F^, the consequences of 



vanishing first variation are more involved (c.f. Lemmas 2.12 - 2.14), and an important 



ingredient is a monotonicity property of double convolutions (c.f. Lemma 2.13| ) which is 
responsible for condition ([3]). Thus, this step holds for any kernel in /C3. (ii) The second 
step involves the explicit parameter regimes and is the reason why we have conveniently 
adopted the hat function. Here we show that there cannot exist a minimiser for Fi, i<2, or 
with fewer interfaces than z if , ([6]) , or ([s]) holds respectively: If a minimiser uq does 
have fewer interfaces than z, then there exists an interval (bounded below in length by w) 
on which z has a bar and uq a space or vice versa. We then contradict our assumption by 
explicitly modifying uq on this interval to achieve lower energy. This last step requires some 
straightforward but tedious calculations. Maple has been a great help in performing the 
many integrations necessary involving the hat function (pa- This step can be reproduced 
for any (pa € IC with different parameter regimes for each specific choice of kernel (pa as a 
result. Specifically, the calculations in the proof of Theorem |1.1| part [2] after (13) or the 
calculations in the proof of Lemma 2.9 after (17) respectively need to be redone for the 
new kernel. 

Collecting the conditions necessary for steps (i) and (ii) we find that it is in principle 



possible to get results of the form of those in Theorem |1.1[ part [2] and Corollary 1.3 
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part [2] for Fi and F2 respectively for any 0o- G /C. Similarly a result for as the one in 
Theorem 1.2 part [2] can be obtained for any (f>a £ ^3- (Note that Corollary 1.3 part [2] 



can be obtained either as a consequence of Theorem |1.2[ part 
latter option allows for less restrictions on (pa- in Corollary 1 1.3 



2] or using Lemma |2.9[ The 



than in Theorem 1.2 



We conclude this section with a few comments on the results, their interpretations and 
limitations. Our brute force arguments are based upon explicit calculations and are as 
such limited to a blurring kernel with small (with respect to uj) support. Indeed, this is 
most discouraging for the deconvolution functionals F2 and -F3 where one one would expect 
the regime of acceptable a to extend far beyond the X-dimension of the underlying bar 
code. Numerical experiments (see Section [s]) support this conjecture. The conditions that 
a, p < w/2 for F2 and F3 may seem particularly alarming since the analogous condition 
for Fi is simply a < oj. However, note that the other condition ^ also puts a restriction 
on the size of a which is essentially of the same form. For F2, one could leave out the 
condition cr < % and not change the principles of the proof, but many more orderings 



in the computation of the integrals become possible (see Remark 2.10) and many more 



calculations need to be done in the proof of Lemma 2.9 Since we still have condition ^ 



on a in place (these extra calculations can only replace ([6]) by a stricter condition on a, if 
anything changes at all) it is doubtful that much can be won by leaving out the condition 



2.14 



cr < ^. For F^ the conditions p, o" < ^ are vital to our proofs via Lemmas 2.12, 2.13 

The numbers in condition ^ may seem a little strange. They are simply a consequence 
of the direct calculations with the hat function. As we have remarked, these calculations 
can be repeated for other kernels in /C. Note that, taking p = co/2, the condition implies 
that the bar code is always recoverable for any a < oj/2, provided A > 20/3uj. 

Another surprising condition might he a < p in Theorem |1.2[ part [2} In general 
if p < (T we do not expect 2; to be a minimiser of F3 over B (or B^^), even if p, o" < 
J and ([5j) are satisfied. A counter example in this case is given hy z = X[o. 425,0.575] 
and u = X[o.425,o.4999] + X[o.500i,o.575] with p = 0.05 and a = 0.06. The fidelity term in 
F3(u) is smaller than the fidelity term in ^3(2;) {\\(pp * u — cp^ * -z||^2(]g) ~ 2.378 • 10~^ 
and \\(pp * z — (pa * -2|li2(]u) ~ 2.407 • 10~^) and thus for A large enough z will not be the 
minimiser of F^ in B. In this particular case the difference is small and so in practical 
applications where A is not too large it might not cause problems, since then the energetic 
cost 2 for two extra interfaces will be much higher than the gain in the fidelity term. Extra 
conditions on the parameters in the case p < a might ensure that u = z is the minimiser 
for F3. The above example suggests that an upper bound for A may be in order. In 
fact, numerical simulations in Section |3] show that the regime p < a poses no problem for 
suitable midrange choices of A. Indeed, they suggest that fixing p comparable with the 
X-dimension and minimising F^ works well for a up to twice the X-dimension. 

We know however that in the degenerate case = p < cr. Theorem |1.1| assures that we 
have u = z as minimiser if both conditions in are satisfied, without an upper bound on 
A. Why the second condition in Q is the correct degenerate form of ^ can be seen by 
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recognizing their common source (16). 

Finally, we note that there is a wealth of work on total variation energy minimisation 
for image analysis. While our results are for rather simple one dimensional images, we 
feel they are novel in that the ID bar code setting entails an image deblurring problem 
of contemporary interest yielding very precise results, and we are unaware of any general 
method for analogous deblurring functionals which would yield similar results. In addition 
to geometric simplicity due to its binary nature (the simplest case of what is called Quantum 
TV in |13|). the bar code problem is different from many other imaging problems in that 
there is a known a priori lower bound on the length scale of the structures in the image 
(via the X-dimension) . An analytically deeper study entails deblurring of 2D bar codes 
121. 



2 Proofs of the Theorems 

2.1 Existence and the trivial minimiser 

Lemma 2.1. Let G IC and fix z ^ (B^ ) and X,a, p > 0. Then minimisers for Fi,F2 
and F3 over B (B^^ ) exist. 

Proof. The proof is a simple application of the direct method in the calculus of variations 
and follows along the same lines for all these functionals. For completeness, we present it for 
Fi and z £ B^j- Let be a minimising sequence for Fi in B, then we can assume every 
Un has bounded L^-norm and bounded BV measure. Therefore, by O §5.2.3 Theorem 4], 
there exists u S BV{[0, 1]) such that Un ^ u in L^{[0, 1]). Since the Un only take values 
and 1 (and a.e. in [0, l]'^), so does u. Thus u e B. 

The total variation is lower semicontinuous under convergence [61 Theorem 1.9] and 
under the special conditions that the functions only take values and 1, so is the 
norm, therefore we conclude via the direct method in the calculus of variations that u is 
a minimiser for Fi. For the functionals F2 and -F3 we use in addition, that the functional 
M I— )• (/>o- * is continuous under convergence, for any a > 0. If we replace B^^ and B by 

and B"^^ respectively the proof does not change. 

□ 

Note that we have not used the fact that (pa is symmetric and unimodal with compact 
support in the above. We only need continuity of u (j)a * u under convergence. 

Next we recall a result about convolutions whose proof follows directly from Fubini's 
Theorem. 

Lemma 2.2. Let f,g,h ^ L^(M) such that f{—x) = f{x), then 

[ [{f*g)-h] = [ [g-{f*h)]. 
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Parts 1 of Theorems 1 1 . 1[ [L2] and Corollary 1.3 follow directly from the following lemma. 

Lemma 2.3. Let i £ {1, 2, 3}, € JC, and A > 0, then the following two statements are 
equivalent: 

1. u = is the unique minimiser of Fi over BV{M). 

2. (a) Ifi = l, IIMU < ^. 

(h) Ifi = 2, Ua*fa\U < h- 
(c) Ifi = 3, Up* fall* < 21- 

Proof The idea of the proof is similar to that in j3 §1.14, Lemma 4] . Note that for general 
u G BV(K) we cannot conclude that (j)^^ * u £ L^(M). For functions u and parameters a or 
p for which this fails we set Fi(u) = oo, i € {2, 3}. 

First let i = 1 . We first prove [T] =^ [2j Assume = is the unique minimiser of Fi in 
BV{R). This is equivalent to, for all h G BV{R) with / 0, 



All/. 



|2 

Il2( 



< 



/ |/i'l+A||/.-/.||i.(i,)= / \h'\ + \mh^^^ + X\\h\\l^^^-2\ [ f^h. (7) 



Because this holds for all h E BV{M.), by rescaling h we can rewrite this as 



2\e / f^h < \e\ / \h'\ + Xe^' 
Jr Jr 



2 



(8) 



for all e G M and all h £ BV{M). Dividing by e, taking the limit e — >• 0, and recognizing 
that e can be positive and negative, we find that (Isl) implies 



fah 



< 



1 

2A 



\h'\, 



for all h£ BV( 



Now per definition 



IIMI* 



sup 

veL^R)J^\v'\<l 



faV 



< 



2A' 



(9) 



where the inequality follows by taking the supremum in ^ over all 
h£{v£ L^{R) : f^v'\ < 1} C BV{R). 

To prove 2 =^ [l let ||/cr||* < j^- Then for all v G L^(M) satisfying \v'\ < 1 we have 



faV 



< 



2A' 



from which it follows that for all h £ BV{i 



IrW\ 



< 



2X' 
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Inequality ^ now follows. 

We proved above that ([T]) implies On the other hand we see that inequality ^ 
implies for h ^ 



[ \h'\+M\f4hm + M\h\\h^^)-2X [ /./i>A||M|i.(^) + A||/i||i,(K)>A||M 



2 



and thus inequality ^ is equivalent to ([T]). This proves the result for z = 1. 

F2 is a special case of -F3 (with p = a). For i = 3 we can derive a statement analogous 
to inequality with h on the left hand side replaced by (f>p * h. Having u = as unique 
minimiser of F3 in BV{M.) is equivalent to 



/ fa ■(l)p*h 



< ^ / \h'\ 
2A 



By Lemma 2.2 we recognise that 



fa- (t>p*h 



R 



for ah h(^BV{l 



* fa -h 



and the result follows as before. 



□ 



Although we assume (pa & '^^ the proof above because that is the most general class 
of kernels we consider, we only use the symmetry and integrability of (pa. 



2.2 Proof of Theorem [TTT], part [2] 



We now turn our attention from the trivial minimiser u = to u = z as minimiser. First, 
we present some elementary results. 



Lemma 2.4. Let (pa^K.,a<h,a<b — a and z = X[a,6]; then 



x£R: fa{x) = H ={a,b} and jx G M : > ^ 



[a,b]. 



Proof. We compute 



fa{a) = (j)a{a-y)dy= (j)a{a-y)dy= (j)^{y)dy=- 

J a J a J —a ^ 

and by symmetry fa{i>) = \- 

Furthermore for x G (a — cr, 6 + a) we compute 

/;(x)= f ct>',{x-y)dy= f (P'a{y)dy+ f c^^dy. 

Ja J x—b Jo 
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The first term on the right is nonnegative if x < b and nonpositive if x > b and the second 
term is nonpositive if x > a and nonnegative if x < a. By symmetry of (pa we then conclude 
that /^(x) > if X < ^ and /^(x) < if x > Moreover we have /^(a) > and 

f'aib) < 0. □ 

Lemma 2.5. Let (pa ^ fC, z ^ Bco, and a < uj, then for every x G M \ suppz we have 

(pa * Z{x) < \. 

Proof. Let x E M \ supp^, then there exist a < b such that b — a > oj, x S {a,b), and 
z{y) = for all y G (a, 6). Define zq := X{-oo,a) + X{b,oo), then 

(pa * z{x) <4>a* zo{x) = f (pa{x - y) dy - I (pa{x - y) dy = 1 - (pa * X[aMi^) < ^• 



The final inequality follows since a < uj < b — a and thus by Lemma 2.4 (pa * X[a,fe] > 5 



(0,6). □ 

The following lemma allows us to consider only minimisers of Fi that have less interfaces 
than z or are equal to z. 

Lemma 2.6. Let (pa £ IC, z €z B^o and let u be a minimiser of Fi over B. Denote by Xi 
the locations of the interfaces of u, then we have for every i 

Consequently if cr <LjJ, Xj is the location of an interface of z for every i. 

Proof. Assume without loss of generality that z ^ 0. Let u minimise Fi over B. We show 
that vanishing first variation ofu implies that at any interface Xj, we must have fa{xi) = \. 
To this end, consider an interface of transition from ?/ = 1 to it = at Xj (the other case 
is treated similarly). By considering a perturbation consisting of extending the u = 1 bar 
up to Xj + i for t small, one obtains no change in the total variation and a change in the 
fidelity term of 

/ [{l-fa?-fl)dx = / {l-2fa)dx. 

Differentiating with respect to t and setting t = gives 1 — 2fa{xi) = 0. 



Let now a < oj. By Lemma 2.4 if z consists of one bar only the i-level set of fa is 



2 

"^^"^ assures 



exactly the set of locations of the interfaces of z. If z has more bars Lemma 2.5 
that the ^-lower level set is not affected by the convolutions of different bars interacting. 

□ 
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Lemma 2.6 tells us that, if a < w, any candidate for minimising Fi over B not equal to z, 
should have less interfaces than z which are located at places where z also has an interface. 
We use this to complete the proof of Theorem First, we introduce a notation that will 
be used frequently in what follows. For a > 0, a < b, and x S M define the functions 

IU.,aM i /' (l ± ^) ,y =!((.- a) (l ± ^) ^ ^) . (10) 

Note that this definition is tailored to the needs of the hat function ^o-. If we want 
to reproduce the calculations that follow for a general convolution kernel G /C we can 
write (j)a- as in ([2| and define /J as 



/•b rb 

I'^{x,a,b) := / p{y — x,a)dy and I^(x,a,a;):= / p{x — y,a)dy. 

J a J a 



We will only state the results for the hat function and hence use the definitions in (10). 



Proof of Theorem 1.1, part [2j Let B 3 uq ^ z he a minimiser of Fi. By Lemma 2.6 



the number of interfaces of uq must be less than the number of interfaces of z and the 
location of every interface of uq coincides with the location of an interface of z. Therefore 
there exists a connected interval C [0, 1] such that |A^| > u; and either 

• z = and uq = 1 on N, or 

• z = 1 and Uq = on N. 

First assume the former case, and let 



u := 



We compute 



Uq on N'^ 

Uq — 1 = z on N. 



u\<2 + / |n'o| (11) 



and 



\U0 - /<7||i2(R) = \\U0 -U + U- fa\\l2f 



\u - + \N\+2 f [uo - n)(u - U), (12) 



where we have used that 

|2 
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Let a £ [0,1] be such that N = [a,a+ \N\] and thus A^^ n [0,1] = [0,a) U (a + |iV|,l], 
then we compute 



(Uo -U){U- fa) 



N 



4>a{x - y)z{y) dydx <2 



a+\N\ 



(j)a{x-y)dy + 



4>a{x - y) dydx 
N yiv=n[o,i] 

(j3tj{x -y)dy \ dx. (13) 



a+\N\ 



Integrals as those in the right hand side of ( 13 ) are commonplace in the proofs of this 
paper. It is therefore very illustrative to work out one of them in detail. Let us consider 

(■a+\N\ f-a 

(paix - y) dy dx. 

Per definition (t)cr{x — y) is zero if |x — y| > a and on its support its value is given by 
^ (1 + iix-y G (-fT,0) and by ^ (l - ifx-y G (0,ct). Let us fix x € [a,a+|7V|] 
for the moment and remember that y G (0, a) in the integral, then (j)a{x — y) 
if 



x-y 



) 



y £ {x — a, x) n (0, a) = < 



ifx-fT<x<0<a, 

{0,x) if X — a < < x < a, 

(0, a) if X — fj < < a < x, 

{x — a,x) if < X — a < X < a, 

{x — cr,a) if < X — a < a < X, 

ifO<a<x-a<x. 



Because x G [a,a + \N\] we can rule out some of these case^and end up with 



y G < 



(0,a) 



if X G (— oo, fj) n (a, oo) 



(x — cr, a) if X G (cr, a + a) Ci (a, oo) 
if X G (a + (T,o+ |iV|). 



1 if o- < a, 

'a, a) if a < a, 
(a, a + a) if a < a, 
(cr, a + cj) if a < cr. 



We see that we have to distinguish between the cases a < a and a < a. Similarly we find 
that 4)a{x — 7/) = ^ (l + ^^) if y G (x, X + fj) n (0, a), which is the empty set because of 
the restrictions on x. 

This now leads us to the computation 



a+\N\ pa 



(j)a{x-y) dydx 



1 

(7 J 
1 



aiA^-'^)dydx 

r/;(i-^) dydx+c'^n^ii 



if CJ < 



dy dx 



if a < 



^For many of the similar calculations in the rest of this paper, a; G R and this kind of simplification will 
not be possible. 
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Because all the integrands are positive, in the case a < o" we can estimate 



1 ) dydx< 



a J x—cr 



la JO 

Therefore we conclude that for both a < a and a < a 

(■a+\N\ f-a 



1 



a 



dy dx. 







<Pa{x — y) dy dx < / IZ{x,x — a,a) dx. 



In a similar fashion we compute 



a+\N\ 



4>aix - y) dydx < 



a+\N\ 



a+\N\-a 



/+(x,a + \N\,x + a) dx. 



During this computation we need to distinguish between the cases a + |A^| < 1 — a and 
a + \N\ > 1 — cj, but as before this distinction doesn't play a role in the final estimate. 



Continuing from (13) we now find 



2 / {uq — u){u — * z) < 2 



IZ {x, x — a,a) dx + 



a+\N\ 



I'L{x,a + \N\,x + a)dx 



a+\N\-(7 



Using this in (11-12) together with |A^| > a; we find 



Fi{u) < Fi(uo) + 2 + Xi^-a-ujj < Fi(no), 

where the second inequality follows by condition ([4]) . This contradicts uq being a minimiser 
of Fi. 

Next we consider the second case, i.e. z = 1 and uq = on A^. We define 



u :- 



uq on N'^ 

uq + 1 = z on N. 



hand side of (12), but for u instead of u: 



As in the first case, we will find an estimate for the integral in the brackets in the right 

\ [ {l-4>,*z) = 2\N\ -2 ! 

JN JN 



JR 



* Z] 



* z. 
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Again we write = [a, a + \N\] and we compute 

(j)„*z= / (f)„{x - y)z{y)dydx > / / (f>„{x - y)z{y) dy dx 
Jn JnJo JnJn 



a+\N\ i-a+lN] ^ 

<Pa{x - y) dydx 



a J a 



i-a+cT i-a+\N\ 

/ IZ{x,a,x) dx + / IZ{x,x — a,x) dx 

J a J a-\-cr 

+ / 11^{x,x,x + a)dx + / I'^{x,x,a+\N\)dx 

Ja Ja+\N\-a 



As in the first case we now find 

Fi{u) < Fiiuo) + 2 + X (^^a - u?j < Fi(no), 
which is again a contradiction with uq being a minimiser. Therefore the only candidate for 



a minimiser is u = z and hence by Lemma 2.1 u = 2; is the unique minimiser. □ 



Remark 2.7. In the above proof everything up to and including (13) is independent of 
the choice of specific blurring kernel and we could have used any (p^^ £ IC. The explicit 
calculations that follow in the remainder of the proof depend on our choice <f>a = (pa, but 
can be redone for a different choice of kernel as explained in the paragraphs preceding the 
proof. 



2.3 Proofs of Theorem [T]2], part [2] and Corollary [Ts], part [2] 

We now turn to F3. 

Lemma 2.8. Let zi,Z2 G B. Both z\ and Z2 on [0,1] consist of a finite collection of 
suhintervals of [0, 1], i.e. alternating bars and spaces. Let ti, i = 1 . . .n and t[, i = 1 . . . n' 
denote the right hand sides of the intervals of z\ and Z2 respectively. In particular tn = 
t'^i = 1. If n > n' , then there exists an interval N C [0, 1] such that [ti, tj+i] C N for some 
i; and for all x G N, either 

zi{x) = and Z2{x) = 1 or zi{x) = 1 and Z2{x) = 0. (14) 

In particular, if zi € Buj, then \N\ >uj. 

Proof. First assume that zi starts with a bar and Z2 starts with a space, i.e. zi = 1 on 
[0,ti] and Z2 = on [0, t'^]. If ti < t[ then [0,ti] C A^. Suppose t'^ < ti. If the conclusion 
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of the lemma is false, then for all i < n' , t[ < ti. This is a contradiction since t^, = = 1- 
If zi starts with a space and Z2 starts with a bar we arrive at a similar conclusion. 

Now assume that zi and Z2 both start with a bar (the situation in which both start 
with a space is similar). Note that zi = 1 on [0, ti] and Z2 = I on [0,t'i]. Suppose ti < t[. 
Then if the conclusion of the lemma is false, we must have t'- < tj+i for i = 1 . . .n' which 
implies 1 = t^, < 1. Suppose ti > t\. If for some i > 1, we have t\ > ti, then the previous 
argument again gives a contradiction. Thus we must have ti > t'- for all i = 2 ... n' — 1. 



But then (14) must hold on one of the intervals for i > n'. □ 



Lemma 2.9. Let z G Bl^ for some i,j S {0, 1}, p^cr < %, (pa = ^o- oind define 



^(-cj3 + 5/9cr2 + 10p3) ifa<p, 
f{p,a):={ ^> PJ J -P. 



Let A, p, and a satisfy in addition 



X 15 

If u £ B^^ is a minimiser of F3 over B^^ , then 



+^[7p + f{p,cj)] <uj. (16) 



\u \ > / \z' 



Proof of Lemma \2.£\ We prove this by contradiction. Let B^^ 3 uq ^ z he a minimiser of 
F3 in B^^ and assume that uq has less interfaces than z, i.e. \uq\ < \z'\. By Lemma 



2.8, there exists a connected interval iV C [0, 1] such that \N\ > u and either 
z = and uq = 1 on A^, or 
z = 1 and uq = on N. 



Define 



then 



and 



u :- 



uq on N'^ 
z on N, 



u\ < / \u'n\ + 2 



*Uo-(f)a* 2:|li2(K) = Up *U-<l)a* Z\\\2(^^) + Up * (mo " m)||^2 

+ 2 / [4>p* {uo - u)) ■ ((f>p*u - 4>a * z 
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from which we conclude that 



F-iiu) <F3{uo) + 2-X{ \\^p*{uo-u)\\l2+2 I Up* {uq - u)) ■ Up* u- c^^ * z) ). (17) 



Because uq — u = ±xn , Lemma A.l gives 



^pHuo-u)\\i2 = \N\--p. 



Next we again distinguish two cases: Case I in which uq = 1 and u = z = on N and Case 
II in which uq = and u = z = 1 on N. We first treat Case I: 



2 / 4>p* {uq - u){(f)p*U- (i>„ * z 



= 2 
> -2 



<pp{x — y)dyl / (f) p{x — w)u{w) dw — j (t>a{x — w)z{w) dw ] dx 

RJN VJiV^nlo,!] JN'^nlO,!] 



(j)p{x — y) dy I (paix — w) dw dx. 

/RJN JAf<=n[0,l] 



Now we subdivide Case I into two subclasses: Case la in which a < p and Case lb in which 
p < a. For Case la we compute 



-2 11 (f)p{x-y)dy[ ^^{x - w) dw dx = 



- bpa'^ - Wp^ 



For details of this computation we refer to (29) in Appendix [A| 
In Case lb the computation is 



R JN 



(f)p{x - y) dy 



(t)a{x — w) dw dx 



15C72 



p^ - 5ap'^ - Wa^ ] , 



the details of which can be found in (30) in Appendix [A| 
In Case II we compute 



p * {uq - U)[(j)p * U - * Z 



= —2 J J (pp{x — y)dy\^ j (pp{x — w)u{w) dw — j (j)a{x — w)z{w) dw^ dx 

> —2 / (l)p{x — y)dyl I (j)p{x — w)dw— I (ft^ix — w) dw ] dx. 
Jr Jn \Jr Jn 

For the first term we find 



-2 [ [ ^p{x -y)dy [ 
Jr Jn Jr 



(u — w) dw dx = —2N. 
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Details of this calculation are given in (31) in Appendix [A} 

For the second term again we need to subdivide into Case Ila in which cr < p and Case 
lib in which p < a. For Case Ila we compute 

(j)p{x -y)dy I 4>a{x - w) dw dx = 2N + (a^ - bpa^ - lOp^j . 



For more details of this computation see (32) in Appendix [A} 

In Case lib we can repeat the calculation with p and a interchanged to get 

1 



2 1 I (j)p{x — y)dy / (p^ix — w) dw dx = 2N + 
IrJn Jn 



15(j2 



p 



5pa^ - Wa' 



Using the combined results of Cases I and II in inequality (17) leads to 

1 



i^3(^) <F3(no) + 2-A \N\ 
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7p + f{p,a] 



<F3(no) + 2 



A 6J 



1 

15 



7/3 + f{p, cr] 



where the final inequality follows from (15) - (16). This contradicts the fact that no is a 
minimiser of F^. 

□ 



2.9 



we have used the conditions z G Bl] , uq £ B^^ , 



Remark 2.10. In the proof of Lemma 
and p,cr < ^ but it might not be immediately clear where. They allow us to order the 
endpoints of the intervals of integration that occur in the integrals in Appendix |A| In 
particular z G Bl^ and no G B^^ imply that the interval N on which z and no differ is 
located at least a distance u away from the endpoints of the interval [0, 1], i.e. a > lo and 
a + I A^l < 1 — uj. If we also take into account the conditions p,cr < ^ we have the ordering, 
for cr < p, 

-p<-a<0<a-p-a<a-a<a<a + a<a + p<a + \N\-p<a + \N\-a 
<a + |iV|<a + |iV|+CT<a + |iV|+p<l-p<l-cj<l<l + o-<l + p 

and an analogous one for p < a. These orderings are important when determining exactly 
which I^{x,a,b) contribute over which x-intervals to integrals like 



K JN 



4>p{x - y) dy 



(paix — w) dw dx. 



'Af<=n[o,i] 

Loosening the condition z £ Bl^ to z £ B^j and consequently no G B^^ to uq £ B is 
possible in principle, but will give rise to more possible orderings of the kind above and 
separate calculations of all the integrals involved need to be done for each possible ordering. 
It is not expected however that this will influence the end result by much if at all. 
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Remark 2.11. Up to and including (17) the steps in the proof of Lemma 2.9 are 
independent of the specific choice of kernels (pa and (pp, but the calculations that make up 
the remainder of the proof do depend on the explicit choice = 4>a- In order to derive 
similar results for other kernels we need to redo those computations with an explicitly given 
alternative choice. 



The result for F2 in Corollary |1.3[ part [2] follows as a direct consequence of Theorem 1.2 



part [2] for by choosing p = a. However, the fact that the fidelity term in F2 vanishes if 
and only ii u = z allows for a direct proof as well. 



Proof of Corollary 1.3, part [2j Since for F2, the fidelity term vanishes at u = z any 



potential competitor must have strictly less interfaces than z. The result follows then 



immediately from Lemma 2.9 with p = a and Lemma 2.1 □ 



To complete the proof of Theorem 1.2 we need a result that tells us that, if cr < /o. 



minimiser of -F3 is either equal to z or has strictly less interfaces. Lemma 2.14 will provide 
exactly this. First we need some preparatory lemmas. 

Lemma 2.12. Let z £ B^, p,a < ^, and ^o- G /C n C(M), then the level-\ set of pp * /o- 
consists of exactly the locations of the interfaces of z. Furthermore the upper level-^ set 
where pp * 4>a * z > ^ is supp z. 

Proof. If z = the results follow trivially. We assume now z ^ 0. First we consider the 
case of a bar code with only one bar. Let a < bhe such that b — a > u and define z := X[a,b]- 
Since the convolution of two symmetric unimodal functions is again a symmetric unimodal 
function (see \14:\ [3] and references therein) we find that (pu* z is a unimodal function with 
mode at xo := ^y^, i.e. (/>o- * z is non-decreasing for x > xq and non-increasing for x < xq, 
and symmetric around x = xq. Since pp is unimodal with mode at x = and symmetric 
around x = we conclude that pp* p^j * z is unimodal with mode at x = xq and symmetric 
around x = xq. Therefore for all x < a and all x > 6 

pp*pa* z{x) <pp*pa* z{a) = pp*p^* z{b) (18) 

and for all x £ [a,b] 

pp*p^* z{x) >4>p*pa* z{a) = cpp*cp^* z{b). (19) 
In the sense of distributions we have 

z =5a- 6b 

where 6x is the Dirac delta measure at x. Hence 

pp*pa* z'{x) = pp* piaix - a) - pip* p^{x - 6). 
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Because (j)p * (pa is unimodal with maximum at we deduce that {4>p -2)' (a) > and 

{(j)p * (pa * zy{b) < 0. Combined with (18) and (19) this imphes that for all x < a and all 
X > b 

(pp* (pa* Z{x) < (pp*(pa* z{a) = (j)p*(pa* z{b) 

and for all x G (a, b) 

(pp*(pa* Z{x) > (pp*(pa* z{a) = (pp*(pa* z{b). 

We now explicitly compute the value cpp* cpa * z{a). 



(pp{z)(pa{q) dq dz 



(Pp * (Pa * z{a) = / / (t)p{a-x)(t)a{x-y)x[a,b]{y)dydx 



a+p rx+a rp 

/ (pp{a - x)(pa{x - y) dy dx = , ^p 

a—p Ja J—pJ—a 
p r-0 r-p r-O 

/ (pp{z)(pa{q)dqdz - / (pp{z)(pa{q)dqdz 

-pJ — a J — pJ—z 

1 I'P 

/ (pp{z)(j)a{q)dqdz. 

p J — z 



In the third equality we have used the change of variables 

-1 



(20) 







+ 



1 -1 



The last equality follows by symmetry of (pa and the fact that <pp and (pia have unit mass. 
Because 



/>0 



(pp{z)(pa{q) dq dz 



pj-z 



we have 



(pp{z)(pa{q) dq dz 



(pp{z)(pa{q) dq dz 



i-Q fp i-z 

I (pp{-z)(l)aiq)dqd{-z) = - 
p Jz Jo Jo 

p rO 

/ (pp{z)(p>a{q)dqdz 

J-z 



fO r-p rO 

(pp{z)(t>a{q)dqdz + / (pp{z)(pa{q)dqdz = {) 



pJ-z 



pJ-z 







and hence by ( 20 ) 



(t)p*(l)a* z{a) 



1 



(21) 



This proves the result if z has only one bar. If z has more bars then we prove that the 
g-lower level set of (pp * fa is the same as the j-lower level set of fa in a similar fashion as 
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the ^-lower level set was identified in the proof of Lemma 2.5 Let x G M \ suppz, then 
there exist c < d such that d — c > uj, x £ {c,d), and z{y) = for all y G (c, d). Define 

zo ■■= X(-oo,c) + X(d,oo), then 

(t>a * fu{x) < (t>p*(l)a* zo{x) = (fip * <f>a{x - y) dy - (j)p* 4>a{x - y) dy 

JR Jc 

= 1 - 0p * 0^ * X[c,d]{^) < 2' 
The last inequality follows from (l)p* (pa * X[c,d] > 5 o^i (c, d) as proven above. 

□ 

Lemma 2.13. Let z £ B^j, 4>p G IC3 and p < ^. Fix x £ suppz, then the function 

{0, p] ^ R : T ^ (/)r * (pp * z{x) 

is non-increasing. 

Proof. First assume that z = X[a,b] for some a < b satisfying b — a> u. We compute 

d f d 

—(j)^*(j)*z{x)=l / (/)p{y-w)—(pr{x-y)dwdy. 

As in ([2]) we write <t)T{x) = p(— t)x[-t,o] (^) + ''")X[o,t] (2^) ^iid thus, by continuity of 
(j)r in T, 

d d d 

:</>r(a;) = X[-Tfi]{x)j7-v{-x,T) + x%T]{x)jrP{x,'r) 



and thus 



d r^'^ d 

— (j)^*cl)p*z{x)= / (j)p{y-w)—p{y-x,T)dwdy 

+ / / (t>p{y-w)-^p{x-y,T)dwdy. 

X — T J Ob 

Using the substitution of variables x'")~'~(o '^l)(u')' 

symmetry of (/^o-, then writing x = x — a and c = b — a and finally dropping the tildes, 
allows us to rewrite the integrals above as the integral in ^ with p instead of a. We can 
thus conclude that 

d 

— (t)r *4'p* Z{x) < 0. 

If z has more bars we note that the double convolution of a single bar of z extends 
a distance oi t + p < 2p < u outside of the bar and thus will not influence the value of 
(j).,- * (j)p * z inside other bars of z. □ 
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Lemma 2.14. Let z G B^^, (po- S IC, o-nd u a minimiser of over B. Denote by Xi the 
locations of the interfaces of u, with xq < xi < . . ., then we have for every i 

4>p * fa{Xi) = ^+4>p*4>p* Ui{Xi), (22) 

where 



Ui :-- 



^ ~ X[xi,Xi+i] if i is even, i.e. Xi is the left interface of a bar of u, 
u — X[xi^i,Xi] if i is odd, i.e. Xi is the right interface of a bar of u. 



(Note that for i even, Ui = tij+ij 

Consequently if (pa = (pa, cr < p < j, and X, p and a satisfy in addition 

^ + (-<y^ + 5pcJ^ + VJp^\ < 00, (23) 



A 15/o2 

then for every i, Xi is the location of an interface of z. 

An example of a bar code u and its accompanying bar codes uq, ui, and U2 is shown 
in Figure [2] 

I I I ■ II 

I I ■ I 

I I ■ 11^ 
I I ■ I 1^ 



Figure 2: For a given bar code u the bar codes uq, ui, and U2 as used in Lemma 
shown 



2.14 



are 



Proof of Lemma 2.14- Let (pa ^K,. Let u minimise F3 over B, then F3 has vanishing first 
variation in u with respect to small perturbations in the locations of the interfaces of u. 
Let xq be the location of an interface of u where the value of u jumps from to 1, in other 
words, it is the left interface of a bar. The argument is analogous for a right interface. We 
consider a perturbed u{t) := uq + X[xo+t,xi\ ^ ^1 where \t\ is small enough such that no 
interfaces are created or annihilated. The number of interfaces of u[t) is equal to that of u 
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and hence we compute (integration is with respect to x) 
\-^[F^{u{t))-F^{u)) 

[(t)p * U{t) - - (^(/>p * -u - 

{4>p * X[xo,xo+t]) - 20p * n • 0p * X[xo,xo+t] + 20p * /a • - u{t)^ 
f {(Pp * X[xo+t,xo]) +'2(t'p*u- (t)p*X[xo+t,xo]+'^(t'p* f<T ■ (u-u{t)^ 



if t > 0, 
if t < 0, 



where we have used that u{t) = u — X[xo,xo+t] if i > and u(t) = u + X[xo+t,xo] if ^ < in 



the last hne as weh as using Lemma 2.2 



Assume for now that t > 0. The case for t < is analogous. Then, using u{t) 

^0 + X[xo+t,xi], 



d_ 
di 



A-i(F3Kt))-F3(n) 



t=o+ 



<t>p{x - y) dy 



d 
) 

' dt 



dt 

xo+t ^ 
(j)p * (f)p * u + 2 



Xo+t 



<t>p{x - y) dy 



XO 



dt 



XO 
XQ+t 



(t^p* fa \ dx 



XO 



= -2(f) p *(f)p* u{xq) - 2(j)p * faixo). 

We can rewrite the first terms as follows: 

</>P * 0p * u{xq) = (f)p*(t)p* {x[xoM + ""ojla^o) = \ + <Pp*(t)p* uo{xo) 



t=o 
(24) 



where we have used (21) to compute (pp* (pp* X[xo,xi]{xo) 



1 



Vanishing of the first variation tells us that the right hand side in ( 24 ) is zero and hence 
1 + 2(t>p* (l)p* Uo{xo) - 2(j)p * fa{xo) = 0, 



which gives equation (22) for Xi = xq. 

Now assume o" < p < ^ and (p^ = (pa- If ^ is such that the white spaces between every 
two subsequent black bars have widths of at least 2p then it follows that (pp*(pp*uo{xi) = 
for every i and equation (22) reduces to 



1 



Lemma 2.12 then completes the argument. Note that in this case condition (23) is not 
necessary. 
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Now assume that u is not as above, i.e. there exist two bars in u separated by a 
white space of width strictly less than 2p. We will show a contradiction. Let xi be the 
right interface of a bar of u and let X2 be the left interface of the next bar, such that 
X2 — xi < 2p < uj. Then the following inequalities should be satisfied 

4>p * fa^Xi) = \ + 4,p* 4>p* Ui{xi) > for i e {1, 2}. 



According to Lemma 2.12 this means that xi,X2 € supp z. Now there are two possibilities. 
The first is that xi and X2 are located in different bars of z. Since z G Bi^ this means that 
X2 — xi > UJ which contradicts our assumption. The second possibility is that xi and X2 
are in the same bar of z. Assume the latter now. 

By the same arguments the right interface of the second bar, i.e. 3:3 also lies in supp z. 
It can lie either in a different bar of z than X2 or in the same one. In the former case we 
have that there exists an interval with |A^| > a; such that z = and n = 1 on A^ and 



using (23) we can use the arguments as in Lemma 2.9 to arrive at a contradiction with the 
fact that n is a minimiser of We conclude that X2 and X3 must lie in the same bar 
of z. In a similar way we find that xq lies in the same bar. If z has more than two bars, 
via induction on the interfaces we find that for every even i, [xi, Xj+i] C suppz. In words, 
every bar of u is contained in a bar of z. 

From the foregoing we deduce that (u — z){x) € {—1,0} a.e. and 

u - z < -X[xi,x2]- (25) 

Define u := u + X[xi,x2]^ then 

o*u- z^ - (j)p*U- 
t>p*u^ + 2(Aa * z ■ (pp* {u- u) - (j)p * uj 

* X[xi,X2]) +2 j^(j>p*U- (i)p* X[xi,X2] - (pa* Z- (f>p* X[xi,X2]) 
(i>p*X[xi,X2]) (j>p*(j>p*U- (f>p*(j)a* Z^ (26) 



where the last equality follows by Lemma 2.2 



■^We don't need u £ S'-* here, because we know that in this construction A'^ has a distance of at least uj 
to a:: = and to x = 1. 



25 



We now use Lemma 2.13 and inequality (25) to estimate 

4'p*4>p*u-4>p*4>„*z] < 4>p* 4>p* {u- z) < - 4'p*4>p*X[xi,x2] 



X2 



Xl 



X2 



0p(x - y)(j)p{y - q) dq dy dx 



Xl 



4>p{y -x)dx\ dy 



<Pp* X[xi,X2] 



Using this in (26) we find 



Fsiu) - F^iu) < -2 - X Up*X[xi,x2]) <0 



which contradicts u being a minimiser. 



□ 



Remark 2.15. The result of Lemma 2.14 doesn't change if z G .B^? and we minimise 
F3 over B^^ for i,j G {0, 1}. Also note that the result can be obtained for any 0o- G /C3 if 



we replace condition (23) by the corresponding parameter range for that choice of kernel, 



which we can obtain be redoing the calculations in the proof of Lemma 2.9 after (17) for 
the new kernel. 



Note that in the case where a < p we could have used Lemma 2.14 in the proof of 
Lemma 12.91 instead of Lemma 12.81 



Proof of Theorem 1.2| , part [2| From Lemma 2.14 it follows that under the stated 
conditions the only possible minimisers of F3 over B^^ are ti = z or a n with strictly less 



interfaces than z. By Lemma 2.9 however such a minimiser cannot have less interfaces 
than z and hence by Lemma 2.1 u = z is the unique minimiser of F^ over B^^ . 

□ 



3 Numerical simulations 

We present a few test simulations for the minimisation problems F2 and F^. To this 
end, there exists an increasing number of state of the art techniques concerning TV-based 
minimisation. However here we are not attempting to write the most efficient algorithm, 
we only aim to test whether the parameter regimes we found theoretically are close to 
optimal or not. Hence we take the naive approach of using a phase field to approximate 
the total variation: That is, choosing e small, we replace the total variation with 
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and consider the LP' gradient descent of the resulting functional. While this technique 
brings in diffuse interface^ (i.e. minimisers will no longer be bar codes), it is well-justified 
for small e (c.f. [I]) in that minimisers will be close to minimisers of the original sharp 
interface problem. One problem with this method in higher dimensions is that one tends 
to get stuck in metastable states, and hence this method would not work well for 2D bar 
codes. However, our ID problem is sufficiently rigid so that the method works well and 
fairly quickly. It takes seconds to run our Python code, and while more direct state of the 
art methods would be substantially faster (as would be needed in a practical application), 
our limited goals are well served by the phase field approach. 
For i<3 the LP gradient flow gives the equation 



ut 

where 



= 2en,, - -W'{u) - 2A 0p * (.^p * n - /,), (27) 



W{u) = " ^^'"^ . In aU of our experiments, a bar code is generated with X- 



dimension u ~ 0.0133. Except for Figure |4] (bottom right), convolution with the hat 
function (pa- is followed by the addition of noise with amplitude a = O.lj^ We used e = 
0.0004 and initial data was always taken to be either u = or u = 1/2. 

The algorithm works well for a far beyond the regime of Theorem |1.2[ We give a few 
sample results. In Figure [3] (left) we see that choosing p = a (i.e. using F2), one obtains 
good results for a larger than twice uj. Figure [3] (right) shows that even for o" ~ 3a;, 
the results are not bad, however they begin to loose accuracy. In Figure |4] we note that 
choosing p to be the X-dimension works well for blurring with a up to twice w. Note that 



here we are in the regime p < a, which we avoided in Theorem 1.2 In fact, our counter 
example suggested that in this regime an upper bound on A is necessary. Figure |4] (bottom 
left) indeed supports this observation by taking A much larger than in Figure [4] (top right). 

We also performed tests where we convolute/blur the bar code with a Gaussian kernel 
with standard deviation a but deconvolute/ dehlur with the hat function (f)p. To obtain 
satisfactory results, one must choose a and p very close to each other and no larger than 
cj, and use a suitably tuned midrange A. We give one example in Figure |4] (bottom right). 

Simulations were also performed for Fi (no deconvolution/deblurring kernel) but we 
always found that using F2 or (with even a small deblurring kernel) was preferable. 



■^For actual implementation, one would need to threshold the output of the minimisation process in order 
to generate a bar code 

'*This choice of constant prefactor in W does not lead to unit surface tension in the sharp interface limit, 
hence A in the simulations differs by an 0(1) multiplicative factor from the A in the analytical results in 
this paper. 

® The added noise was determined as follows: The 400 grid points making up each interval of length uj 
were divided into 16 equal groups, each of which was assigned a random number between —a and a. 
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Iambda=100000, sigma=0.0300, a=0.10 
omega=0.0133, rho=0.0300, epsilon = 0.0004 





lambda=100000, sigma=0.0400, a=0.10 
omega=0.0133, rho=0.0400, epsilon = 0.0004 







Figure 3: Here we look at minimisers of F2 {p = a) to find that the algorithm works for 
blurring far past w. In all simulations, the three rows are as follows: A bar code is generated 
with X-dimension exactly a; = 0.0133; convolution f„ of the bar code with 4>a- with added 



noise of amplitude a; final steady state for (27) superimposed with the generating bar code 



4 Discussion 



We have presented results on the accuracy of TV-based energy minimisation methods for 
bar code deblurring in certain parameter regimes. Numerical simulations, which included 
the effects of noise, show that these methods are valid in much larger regimes and in par- 
ticular, allow for significantly more blurring. While our analytical results did not showcase 
the benefits of using a deconvolution/deblurring kernel (i.e. the merits of F2tF^ versus 
Fi), numerical experiments showed clearly that the presence of a deblurring kernel in F2 
or F3 always gave better results over no deblurring {Fi). 

In practice, the size of the blurring kernel pertains to the so-called spot diameter of the 
laser beam at impact with the bar code. This spot diameter is a function of the laser beam 
and the distance from the scanner to the bar code. According to Palmer |10j (p. 127), most 
scanners can successfully read a bar code if the spot diameter is no greater than \/2 times 
the X— dimension (i.e. for 2a < \/^uj). This suggests that the even the conditions we have 
imposed on a in our results are not completely unreasonable. However, as suggested by 
the numerics, one might be able to prove results for a past the X— dimension. 

Realistically neither the spot diameter (size of the blurring kernel) nor the distribu- 
tion of the beam intensity (shape of the blurring kernel) is exactly known, and inferring 
this information from signals is an ill-posed problem. In [3], the author considers a Gaus- 
sian ansatz for all kernels but introduces a novel optimization scheme for determining the 
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Iambda=1000, sigma=0.0220, a=0.10 
omega=0.0133, rho=0.0133, epsilon = 0.0004 



lambda=1000, sigma=0.0260, a=0.10 
omega=0.0133, rho=0.0133, epsilon = 0.0004 





0.2 0.4 



lannbda=10000, sigma=0.022O, a=0.10 
omega=0.0133, rho=0.0133, epsilon=0.0004 



1.2 
1.0 
0.8 
06 
0.4 
02 
oo I — I 




convolution with Gaussian, lambda=800, sigma = 0.0130, a=0.10 
omega=0.0133, rho=0.0130, epsilon = 0.0004 

1.2 F ' ' • • 




Figure 4: Minimisers for F3. Here we take the deconvolution kernel size to be the X- 
dimension. Top row: we convolute the data with X-dimension exactly lo for two choices 
of a. The algorithm works well with A = 1000. However, as noted in the bottom left, 
for larger values of A it loses information. Bottom right: a bar code is convoluted with a 
Gaussian with standard deviation a but deconvoluted with 6a ■ 
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standard deviation of the blurring kernel. In terms of the shape of the kernel, our last 
simulation in Figure |4] (bottom right) is suggestive. We note that if the convolution in the 
measured signal is done with an infinitely supported Gaussian with standard deviation a, 
then deconvolution with a hat function of approximate size a works reasonably well. Thus 
if one could determine certain statistics of the blurring kernel, one could then deconvolute 
with a set kernel possessing similar statistics. Determining such statistics should in princi- 
ple be possible as some standard bar code symbologies have a fixed structure at their left 
and right boundaries (the left and right guards, c.f. [lOj). 
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A Calculations in the proof of Lemma 2.9 



In this appendix we collect some of the longer calculations in the proof of Lemma 2.9 We 
start with a lemma. 

Lemma A.l. Let z := X[a,h\ for some a < b, a < and /o- = <^o- * z, then 

fl = b-a-^-a. 
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Proof. We compute 

r 



if X G (—00, a — a], 

I^{x,a,x + a) if X G [a — fj, a], 

IZ{x, a, x) + /^(x, x,x + a) if x G [a, a + a], 

(x, X — cj, x) + (x, X, X + cj) if X G [a + a,b — a], 

IZ{x, X — a, x) + I!^{x, x,b) if x G [6 — cr, 6] , 

IZ{x,x-a,b) if X G [6, fe + cr], 

if X G [6 + (7, 00). 



r 



K{x + (T-a) 



2(72 



if X G (—00, a — a], 
if X G [a — fj, a] , 



^[x-a-{l + V2)a)(x-a-{l-V2)cr) if X G [a, a + o"] 



1 



if X G [a + a,b-a], (28) 



1 



2(t2 

^{x-a-by 




X - 6 + (1 - \/2)cTj (^x - 6 + (1 + ^/2)cjj if X G [6 - a, 6], 

^2 [ixe[b,b + a], 

if X G [6 + o", 00). 

An explicit computation of the integral we are interested in leads to the result. 



□ 



Next we give the calculations for the different cases described in the proof of Lemma 2.9 
Case la: 

(j)a{x — w) dw dx 



2 / / 4)p{x -y)dy [ 



/^(x, a, x + p) I (x, X — cr, x) + /+(x, x,x + a) \ dx 



+ J /^(x, a, X + p) ^/^ (x, X — (T, x) + /^(x, X, a)^ dx 
+ J ^I'^{x,a,x) + Il{x,x,x + p)^IZ{x,x — a,a) dx 

fa+\N\ / \ 

+ / /^^(x,x-/9,x) + /:^(x,x,a + |iV|) /+(x,a + |7V|,x + cr)(ix 

Ja+\N\-a\ J 
f-a+\N\+a / \ 

+ / I^{x,x - p,a + \N\)[ r{x,a + \N\,x) + I'^{x,x,x + a)] dx 

Ja+\N\ \ J 



+ 



ra+\N\+p / \ ^ 

/ I'L{x,x- p,a + \N\) IZ{x,x- a,x) + /+(x,x,x + a)) dx> 

Ja+\N\+a V / J 



- 5pa^ - Wp^ ] . 



(29) 
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The way to find the specific intervals of integration in the integrals above (and those 
that follow below) is similar in spirit to what was done in the proof of Theorem We 
do not give all the details here, but it is important to reflect on the role of the conditions 
z G i3^, no G (instead oi z ^ and u G B) and p,(7 < ^. Such considerations are 
addressed in Remark |2.10[ 

Case lb: 



-2 / / (f)p{x — y) dy / (paix — w) dw dx 

Jr Jn JAf=n[o,i] 

-2|y I+{x, a,x + p) ^I'^{x, X — a,x) + I^{x, x, a)^ dx 

+ J ^I^{x,a,x) + I'^{x,x,x + p)^IZ{x,x — a,a) dx 



+ / [ I^{x,x — p,x) + I'l{x,x,x + p) ]lZ{x,x — a,a) dx 

J a+p 
fa+\N\-p 



ra+\N\-p / \ 

+ / [l'^{x,x — p,x) + l'^{x,x,x + p)\l'^{x,a+\N\,x + a)dx 

Ja+\N\-a \ J 
ra+\N\ / \ 

+ / [l''_{x,x - p,x) + l1{x,x,a + \N\)\l1{x,a + \N\,x + a)dx 

Ja+\N\-p\ J 
l-a+\N\+p I \ 

+ / lt{x,x- p,a^\N\)[r_{x,a^\N\,x)^l%{x,x,x^a)\dx 

Ja+\N\ \ J 



'a+|Af| 

^ ' p^ - bap^ - lOo-^ ) . (30) 



15^2 
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Case II, first term: 

— 2 / / <i)p{x — y)dy / (j)p{u — w) dw dx 
Jr Jn Jr 

= J -^+(^' 3; + p) ^I'!_{x, X — p, x) + I^{x, x,x + p)^ dx 

J ^I'L{x, a, x) + I'^{x, x,x + p) dy^ ■ ^I^{x, x — p,x) + I^{x, x,x + p)^ dx 



+ 



+ / I X — p, x) + /^(x, X, X + p) I • I /^^(x, x — p, x) + I^!.(a;, X, X + p) J 

+ / f/^^(x, X — p,x) + /^(x, X, a + |Ar|)V ('/^(x, X — p,x) + /^(x, X, X + p)^ dx 

Ja+\N\-p\ J \ J 

/■a+|Af|+p / \ ] 

y 7^(x,x-p,a + |iV|)|^7^(x,x-p,x) + /^(x,x,x + p)J dx\ 



+ 

= -2iV. (31) 
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Case II, second term, Ila: 



2 / / (ppix — y) dy / (paix — w) dw dx 
Js. Jn Jn 



= 2< / I^{x,a,x + p)I'^{x,a,x + a) dx 

I J a— a 

+ J ^I^{x,a,x) + I'^{x,x,x + p)^ ■ ^IZ{x,a,x) + I'^{x,x,x + a)^ dx 
+ J ^l'!_{x,a,x) + I^{x,x,x + p)^ ■ ^IZixjX — a,x) + I'^{x,x,x + a)^ dx 



+ / { I^{x,x — p,x) + I^x, x,x + p) \ ■ I I'L{x,x — a,x) + I^{x, x, x + a) \ dx 

fa+\N\~a / \ / \ 

+ / { I'L{x,x - p,x) + Il{x,x,a+\N\)) ■ [ I'L{x,x - a,x) + I'^{x,x,x + a)] dx 

Ja+\N\-p V /V / 

+ / ' ' ( I^{x,x - p,x) + Il{x,x,a + \N\)] ■ ( IZ{x,x - cr,x) + Il{x,x,a + \N\)] dx 

Ja+\N\-a\ J \ J 

l-a+\N\+a 

+ / I^{x,x- p,a+\N\)r{x,x-a,a + \N\)dx 

Ja+\N\ 



= 2N + ^(^a^ -5pa^ -Wp^y (32) 

B Proof that (j)^- satisfies condition (3) 

In this Appendix we prove that the hat function cj}a- satisfies condition We do this by 
proving a more general result first and then showing that this holds for the hat function 
in particular. 

We use the notation as in ^ and introduce 

Lemma B.l. Use the notation as in If for each r G (0, cr] 
d 

1. either ^p(x, r) is monotonically increasing in x and J ((J., r, 0, c) < for all c > 2a, 

OT 



d 

2. or -^p(x,t) is monotonically decreasing in x and J'(a,T, §,c) < for all c > 2a, 

OT 



then J'{a,T,x,c) < for all r G (0,o"], for all c > 2a and all x G [0,c], i.e condition ^ 
holds. 
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Proof. Let o" > 0, c > 2a, and r G (0, a]. 
Define fa := (^<t * X[o,c] and 




if - r < X < 0, 
if < X < T, 
otherwise. 



This allows us to rewrite 



We first consider case 



J{a,T,X,C) =1pr* fcr{x). 

d 



1 Since -g-p{-,T) is monotonically increasing, the function — V't 

is symmetric unimodal. "Since the convolution of two symmetric unimodal functions is 
again a symmetric unimodal function (see [14^ [3] and references therein) we find that 
— V'r * /o- is a unimodal function with mode at |. Hence ipr * /o-(O) = ipr * /cr(c) and 
A * fa{x) < Vr * fa{0) < for all X G [0, c]. 

Next we consider case [2} In this case ipr is symmetric unimodal and hence ipr * fa is 
unimodal with mode at | and hence for all x G [0, c] we have ipr * faix) < ipr * faic/2) < 
0. □ 



We complete the proof that (pa- satisfies condition ([3|) by showing that (pa satisfies con- 
dition[l]in Lemma B.l Let a, r, c and x satisfy the conditions in ([sj). For the hat function 

1 / x\ d 1 / 2x \ 
we have p{x,t) = — ^1 j and hence —p{x,T) = — ( —1 H ) is monotonically in- 
creasing in X. Furthermore by (28) — with a 
have fa{y) = 1 — fa{—y)- We then compute 



\ T 

0,6 = c — we find that for y G [0, r] we 



J(cT,r,0,c) = Vr*/a(0) 



-1 + 



1 

22/ 



-1 + 



2y 



{fa{-y) + 1 - fa{-y)) dy 



T 



dy = 0. 
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